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t j ! Abstract: We investigate the response function of hot nuclear matter to a small isovector 

^ [ external field using a simplified Skyrme interaction reproducing the value of the symmetry 

qJq ■ energy coefficient. We consider values of the momentum transfer corresponding to the dipole 

§ . oscillation in heavy nuclei. We find that while at zero temperature the particle hole interaction 

O ' is almost repulsive enough to have a sharp (zero sound type) collective oscillation, such is no 

On ' longer the case at temperatures of a few MeV. As a result a broadening of the dipole resonance 

occurs, leading to its quasi disappearence by the time the temperature reaches 5 MeV. The 

O ' sensivity of the temperature evolution of the width when modifying the residual interaction 

^ ' strength is also examined. 
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A large amount of data has been accumulated over the years concerning the properties 
of giant resonances in hot nuclei |jl], and the evolution of their energy and width as tem- 
perature increases. To a good accuracy the energy of giant dipole resonance is found to be 
independent of the excitation energy while its width is found first to increase at low excitations 
energies with a saturation at higher energies 0. On the theoretical side the energy of the 
giant dipole resonance is well understood already at the level of the random phase approx- 
imation (RPA). Concerning the increase of the width with T however, the situation is still 
unsettled. Indeed although there are various mechanisms contributing to this increase such as 
shape fluctuations or angular momentum effects |^ their total contribution does not seem 
to account for the observed data. Furthermore there are also dynamical effects inhibiting this 
increase . It is worthwhile noting that in some of the early RPA calculations using oscillator 
functions |0 a rapid increase of the width was found. This result was however not confirmed 
in the self-consistent RPA calculations including continuum states performed by Sagawa and 
Bertsch ^j. Recently the early RPA calculations have been reexamined by Nicole Vinh Mau 
who concluded that one important ingredient in these calculations was the appearence of 
new particle-hole configurations as temperature rises, an effect to which more attention should 
be given when trying to explain the broadening of the giant dipole resonance. Why are these 
effects found to be nearly negligible in the calculations of Sagawa and Bertsch? The aim of 
the present letter is to examine this question in the framework of a self consistent calculation 
of the response function of nuclear matter using a schematic Skyrme force. The advantage of 
this calculation is that while incorporating a full self-consistency it is still tractable and leads 
to transparent formulae. In particular these formulae show that while at zero temperature the 
isovector particle -hole residual interaction is almost repulsive enough to have a well developped 
zero sound, such is no longer the case at temperatures of a few MeV, leading in general to a 
rapid weakenig of collective effects. 

In the case of nuclear matter single particle states are labeled by their wave number k, their 
spin 0" and their isospin r. We wish to calculate the response to an external field of the form: 

Kxt = er3e^*'i-''e-^(-+^'')*. (1) 



2 



where Ts is the third component of the isospin operator and t] a vanishingly small positive 
number corresponding to an adiabatic switching of the external field. For this purpose we use 
mean field theory. In this case the evolution of the one-body density matrix p is determined by 
the time dependent Hartree-Fock equation: 

tdtp=[W + Ve.t,p], (2) 

where W is the mean field hamiltonian. Let us consider for simplicity the case of a schematic 
Skyrme force: 

V = to(l + XoPa)S{ri - ra) + h6{ri - r2)5(r2 - rg). (3) 

In this case we have 

= |1 + U(r), (4) 
where [/ is a local potential, given in the case of neutrons by: 

t/„(r, t) = |{2p(r, t) - p„(r, t)} + |{p2(r, t) - p^r, t)}, (5) 

with a similar expression for protons. In the previous equation p(r, t) is the density distribution 
< r|p|r >. An identical equation would be obtained with a density dependent version of the 



Skyrme force |]T0| and thesubsequent developpement would be unmodified. For a small enough 
external field it is legitimate to linearize the mean field evolution equation (^. This procedure 
leads to the following approximate equation for the difference 6p = p„ — Pp of the neutron and 
proton density matrices: 

tdt < k\5p\k' >= {e(k) - e(k')} < k|5p|k' > +{/(kO - /(k)} < k|(f/„ - Up)\k' > 

(6) 

+2e{/(k') - /(k)}5(k' - k - q)e-^('^+^'7)*. 

In this equation e(k) = h'^k^/2m is the energy of the single particle state with wave number k 
and 

/(k) = l/{l + e^(^«~^)}, (7) 

is the corresponding occupation number. Note that in the linear response approximation we 
have from equation (H): 

U^ir,t)-Upir,t) = 2Vo6pir,t) (8) 



where 



^o = -|(xo + ^)-|po, (9) 



Po being the saturation density of nuclear matter. Note that Vq is related to the symmetry 
energy coefficient a^- of nuclear matter via the relation: 

ar = ^Tf + Wopo. (10) 

where Tp is the kinetic energy at the Fermi surface. The above equations suggest to look for a 
solution of the evolution equation (^) in the form: 

< k|(5p(t)|k' >= (5(k' - k - q) < k\5pit = 0)|k + q > e-'("+*'')*. (11) 

For this Ansatz the change in the density distribution has the following structure: 

5p(r,t) = ae-'^-'"e-'("+''')*, (12) 

where: 

a = ^ < k|5p(t = 0)|k + q > . (13) 

k 

The change between the neutron and proton potentials is thus found from equation (|^) to be: 

< k\{Un - Up)\k' >= 2Voa6{k' - k + q)e-^('^+*'^)*. (14) 

Inserting this expression into the linearized evolution equation for the density matrix (|^) we 
find that our Ansatz is indeed a solution provided that: 

' ^ ^' LU + IT] - e{k) + €{k + q) ^ ' 



Returning to equation (|13]) we obtain a linear relation determining the value of a: 

a = {Voa + e)IlQR{uj, q), (16) 

In this equation Iloij is the unperturbed retarded response function defined by: 

n r ^ 2 f^^^ /(k + q)-/(k) 

noH^,q = jTT^ / dk — — — -, (17) 

(27r)'^ J uj + 17] — e(k) + e(k + q) 
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Solving equation (|TB|) for a we obtain the following expression for the retarded response function 



Ilij in the RPA approximation: 

e 1 - VonoHlc^, q) 

The result we have derived is valid only in the case of a simplified Skyrme force. For other 
interactions the quantity Vq would not be a c-number. The case of a full Skyrme force was 



worked out at zero temperature by Garcia- Recio et al |jTl| who still obtained a closed expression 
for the response function at the cost of introducing additional unperturbed functions. For a 
given momentum q we find that there is a resonant response when the frequency lo corresponds 
to a zero of the denominator, i.e., when: 

l = yonoij(u;,q). (19) 

The real part of lo determines the energy of the collective mode while '^m[ijj) determines its 
life time . The importance of the resonance condition (^) is also seen when looking at the 
distribution of strength per unit volume S[iS) for the operator exp(iq.r). Indeed it is given by 



-S^mn/j/vr i.e. [p 



For symmetric nuclear matter this function enjoys the following sum rule 

/•CO fp' 

ujS{uj)duj = ^Pol^- (21) 



Let us now try to apply the previous formulae to the case of finite nuclei. For this purpose 
a usefuU guide is the Steinwedel and Jenssen model [|1^]. In this model neutrons and protons 
oscillate inside a sphere of radius R according to the formula: 

p„(r, t) - pp{r, t) = e sin(q.r) exp{iujt), (22) 

the total density remaining equal to the saturation density po and the wavenumber q being 
given by 
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In order to describe the dipole resonance in finite nuclei this model thus suggests to look at the 
strength function S{uj, q) calculated at a momentum transfer q = 7t/2R, where R is the nuclear 
radius. In the following we shall focus on the nucleus lead-208 for which we take R = 6.7 fm 
and q = .23/m~^. We use the following values of the parameters: to = —983.4 MeV x fm^, 
^3 = 13106 MeV X fm^ and Xq = .48. These values were fitted to reproduce the binding energy 
(E/A= -16MeV), the saturation density( po = 0-17 fm"'^) and the symmetry energy(aT- = 30 
MeV) of nuclear matter. For these values Vq = 203MeV fm^. 

To calculate S{uj) we still need the expression of the unperturbed Lindhard function. The 
imaginary part can be expressed in terms of elementary functions as |T^: 



inq 

where (3 = 1/ kT while 



^ ^ , , kT f2m\\ l + e^(^+"/2) 

^mUoRiu,q) = -—^j^) \og^-^-^^^-^ (24) 



A = f.-^-'M (25) 
^ 4e(g) 4 ^ ^ 

In this equation fi is the chemical potential and e(g) the single particle energy of the state with 
momentum q. The real part can be expressed as: 

^eUoR = -J F{k)df{k) (26) 



where f{k) is the occupation number while: 

, 1 mk \ ^ k 



,mw ^ q ^ ^ ' 

'hkq 2k hkq 2k' 



(27) 



with: 



0(a;) = (1 - a;)(l + x)log|^^|. (28) 

Note that at zero temperature we have df{k) = —6{k — kF)dk so that F{k) is just the zero 
temperature value of the real part of Hqr. This implies that equation (B^) is well suited for a 



numerical integration centered around the value k = \j2m^/h^. 

The real part of the Lindhard function is graphed on figure 1 as a function of the frequency 
uj for q = .23/m~^ and several values of the temperature T=0,2,4,6 MeV. It can be seen that 
at T=0 the quantity ^eUoR has a maximum near u = ISMeV which is not quite sufficient to 
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produce a zero in the quantity 1 — VoS^elloR. Such a zero would correspond to a sharp collective 



oscillation analogous to the zero sound mode in neutral Fermi liquid systems |12|. Still, for 
T=0, the value of VoSJello/j is sufficiently close to unity at its maximum to produce a peak in 
the strength at u = lAMeV, as seen in figure 2. It should be noted that the position of this 
peak is in close agreement with the observed value (14 MeV). It is also interesting to note in 
figure 1 that as temperature increases the maximum in ^cUqr becomes weaker and weaker as a 
result of the averaging with the occupation numbers in equation 26. This leads to a broadening 
of the peak in the response function as can be checked in figure 2. By the time one reaches a 
temperature of 6 MeV only small collective effects remain visible. 

The analysis we have just given also explains why some RPA calculations do give a rapid 
increase of the dipole width and some others do not. Indeed a slightly stronger residual interac- 
tion would be able to produce at zero temperature a sharp zero sound type mode which would 
disappear rapidly as temperature increases. In contrast a weaker residual interaction gives at 
zero temperature an already broad peak evolving slowly with temperature. An illustration of 
this point is given in the lower part of figure 2, which shows the evolution of the dipole mode for 
a residual interaction with Vq = 100 MeV x fm^. Similar results are obtained for lighter nuclei 
i.e. larger values of the momentum transfer q. One difference however is that the maximum of 
the real part of the response function becomes less pronounced as q increases which produces 
broader resonances. Of course the discussion we have given ignores finite size or surface effects 
and is further limited to the RPA framework. Some effects not included in this framework, 
such as the coupling to two particle-two hole states, are known to be important to describe 
the damping of collective vibrations |T^. Our analysis was further based on the Steinwedel 
and Jenssen model. We believe that analogous results would be obtained with the Goldhaber 
Teller model Indeed we expect in this case the response function to be dominated by the 
q = Tx/2R mode with small corrections only arising from the multiples q = nir /2R {n > 1), 
because the peak in ^eUoR is much smaller for n > 1. 
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Figure Captions 



Figure 1 Real part of the non perturbed Lindhard function multiplied by he (in fm ^) as a 
function of the energy uj (in MeV) for a momentum q=.23 fm~^ and for different values of the 
temperature T=0,2,4 and 6 MeV. 

Figure 2 Distribution of strength per unit volume for the operator exp(iq.r) (in MeV~^x 
fm~^) as a function of the energy u (in MeV) for a momentum q=.23 fm~^ and for different 

values of the temperature T=0,2,4 and 6 MeV. The upper curves correspond to an interaction 
strength Vo= 203 MeVxfm^ while the lower curve corresponds to Vq— 100 MeV x fm"^. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9404008v2 
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